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Abstract 

We study the finite quark mass effects of the low-energy QCD using the improved lad- 
der Schwinger-Dyson and Bethe-Salpeter equations which are derived in the manner con- 
sistent with the vector and axial- vector Ward-Takahashi identities. The non-perturbative 
mass-independent renormalization allows us to calculate the quark condensate for a non- 
zero quark mass. We explicitly show that the PCAC relation holds. The key ingredients 
are the Cornwall-Jackiw-Tomboulis effective action, the generalized Noether current and 
the introduction of the regularization function to the Lagrangian. The reasonable values 
of the pion mass, the pion decay constant and the quark condensate are obtained with a 
rather large Aqcd- The pion mass square and the pion decay constant are almost propor- 
tional to the current quark mass up to the strange quark mass region. It suggests that 
the chiral perturbation is applicable up to the strange quark mass region. We study the 
validity of the approximation often used in solving the Bethe-Salpeter equations too. 
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1 Introduction 



The program to derive the observed properties of hadrons non-perturbatively in QCD has 
been pursued with great intensity but not accomplished yet. The concept of chiral symmetry 
and its spontaneous breakdown are among the most important aspects of low-energy hadron 
physics. The spontaneous breakdown of chiral symmetry is believed to be responsible for a 
large part of the low-lying hadron masses as well as for the emergence of octet pseudoscalar 
mesons as Goldstone bosons. In order to explain the observed hadron spectrum, one also needs 
small, explicitly chiral-symmetry breaking terms, namely, the flavor-dependent current quark 
mass terms. 

The Conwall-Jackiw-Tomboulis (CJT) effective action approach ]]J for composite opera- 
tors is widely used to study the dynamical symmetry breaking phenomena in the quantum field 
theories. The extremum condition for the effective action with respect to the quark propagator 
leads to the Schwinger-Dyson (SD) equation for the quark propagator on the non-perturbative 
vacuum and the second variational derivative of the effective action with respect to the quark 
propagator leads to the Bethe-Salpeter (BS) equation describing the bound states. The ad- 
vantage of the present approach is that the derived SD and BS equations are consistent with 
the symmetry of the effective action evaluated in a certain approximation scheme. 

The QCD SD equation for the quark propagator has been studied in the improved ladder 
approximation (ILA) by Higashijima [fj| and Miransky [[|. They took the ladder diagrams 
of one-gluon exchange between q and q and assumed that the coupling constant is modified 
according to the standard perturbative corrections. It has been shown that the asymptotic 
behavior of the solution is consistent with the leading order renormalization group analysis 
while the infrared gluon exchange breaks chiral symmetry dynamically. Aoki et al. solved the 
BS equation the J PC = 0~ + qq state and confirmed the existence of the Nambu-Goldstone 
pion in this approximation M. The numerical predictions of the pion decay constant f n and 
the quark condensate (ipip) are rather good. It was also shown that the BS amplitude shows 
the correct asymptotic behavior as predicted by the OPE in QCD [||. The masses and decay 
constants for the lowest lying scalar, vector and axial-vector mesons have been evaluated by 
calculating the two point correlation functions for the composite operators ipMip. The obtained 
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values are in good agreement with the observed ones ||. 

So far, the current quark mass term has not been introduced in the studies of the BS 
amplitudes in the ILA. The purpose of the present paper is to investigate the effects of the 
finite current quark masses on the BS amplitudes for the J PC = h states. As shown in 
H, the asymptotic behavior of the solution of the SD equation for the quark propagator with 
the finite current quark mass is rather different from that in the chiral limit. Therefore, it is 
important to study the effects of the finite current quark masses not only on the SD equations 
but also on the BS amplitudes. 

There have been many studies of the pion BS amplitude using the effective models of QCD 
and/or the approximation schemes of the QCD |7], RJ. The advantages of the ILA model are 
as follows, (i) The model is given in the Lagrangian form so that one is able to apply the CJT 
effective action formulation and study symmetry properties of the system, (ii) The asymptotic 
behavior of the solutions of the SD and BS equations is consistent with the renormalization 
group analysis of QCD. (iii) It has been shown that the ILA model corresponds to the local 
potential approximation with the ladder part in the non-perturbative renormalization group 
approach ||. (iv) The angular integration in the SD equation can be performed analytically. 
On the other hand, the disadvantages of the ILA model are as follows, (i) The axial-vector 
Ward-Takahashi identity is violated [ iPf] . (ii) The quark may not be confined in the color 
singlet state. 

In the finite quark mass case, it has been known that there is a difficulty in defining the 
quark condensate in the studies of the QCD SD equation. The extraction of the perturbative 
quark mass contribution in the UV region is not sufficient to remove the UV divergence in 
the SD equation since the SD equation includes the non-perturbative contribution even in the 
UV region. Inspired by Kusaka's idea of the non-perturbative renormalization of the fermion 



mass term in the mass-independent renormalization scheme pT| , we propose a novel way to 
renormalize the SD equation and define the quark condensate with the finite quark mass. 

The paper is organized as follows. In Sec. |2| we explain the ILA model Lagrangian we have 
used in the present study. In Sec. |3] the SD equation is derived from the CJT action and the 
renormalization of the SD equation is discussed in Sec. [|. In Sec. [5] our method of solving 
the BS equation for the pseudoscalar meson is presented. In Sec. |6| the formulation for the 
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meson decay constant is given and the low-energy theorem is discussed. Sec. [F] is devoted to 
the numerical results. Finally, summary and concluding remarks are given in Sec. 



2 Improved Ladder Model of QCD 

We work with the following Lagrangian density of the improved ladder approximation(ILA) 
model of QCD proposed by Aoki et al. j|, |6|, 

C[tp,tp} := C iTee [ilj,i/j] + C int [ip,ip} , (1) 

CfreeWM ■= W ~ ™o)V> ■ (2) 

Here the function f(() of ( = d 2 is introduced to provide a cut-off regularization of the ul- 
traviolet divergences of the quark loops. The reason we introduce the cut-off function at the 
Lagrangian level is to preserve the consistency between the SD and BS equations. If one uses 
the regularization that is inconsistent between the SD and BS equations, the low-energy rela- 
tion based on the chiral symmetry should be violated by the regularization. The function f(() 
should satisfy f(( = 0) = 1 and /(£) — > oo for ( ^> A^y In this paper, we employ the sharp 
cut-off function 

/(C) = l + M0(C-A*v), M^oo. (3) 

We introduce the bare mass of quarks mo which is evaluated at Auv-[f|) ^ n g enera l m o is a 
diagonal flavor matrix i.e. mo = diag(m u , rrid, m s ) for Nf = 3. In this paper we deal only with 
a flavor independent mass and therefore the case with SU(3)f symmetry. 
The interaction term is given by 

AnthMK*) := -\l JC mm '> nn '(p,p';q,q>) 

Z Jpp'qq' 

X VM^WAMMfl*^*^ , (4) 

xW ^ (l±Pl. _ i±A ^^T a ) mm \ lu T a ) nn ' (5) 

.4 

where f p denotes / an d Pe represents the Euclidean momentum. The Fourier transfor- 
mations of fields are defined by ip(p) — J d 4 xe ipx ip(x) and ip(p) — J d 4 xe ipx ijj(x). The indices 



m, n, • • • are combined indices m := (a,i, f), n := (b,j, g), ■ ■ ■ with Dirac indices a,b, ■ ■ ■ and 
color indices ■ ■ ■ and flavor indices f,g,- • ■■ T a denotes the generator of the color SU(N C ). 
According to Higashijima and Miransky, we choose a particular set of the momenta that deter- 
mines the running coupling constant, i.e., 



?{p E , ?e) = We - 4) Ape) + 0(q E - PeWe) 



(6) 



This way of introducing the running coupling constant is very natural from the non-perturbative 
renormalization group approach with the local potential approximation |J. It is often called 
the Higashijima- Miransky approximation. The infrared cut-off tjp is introduced in the running 
coupling constant as 

1 1 
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IF 



(7) 



t := In 



Po- 



Pe 



A 2 



1 3tiF — to -\- 2 
{ 2fo (1+tiF) 2 

1, 



for t < t 



1 11N C - 2N 



(9) 



(4tt) 2 3 

Above t IF , g 2 (p%) develops according to the one-loop result of the QCD renormalization group 
equation and below to, 9 2 (Pe) * s kept constant. These two regions are connected by the 
quadratic polynomial so that g 2 {p 2 E ) becomes a smooth function. Here Nq is the number 
of colors and Nf is the number of active flavors. We use Nc = Nf = 3 in our numerical studies. 
The gluon propagator is given in the Landau gauge 



iD^ik) = [g^ - 



k"k l 



A- 2 



(10) 
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3 SD equation 

In order to derive the Schwinger-Dyson (SD) equation, we use the formalism of the Cornwall- 
Jackiw-Tomboulis (CJT) effective action which is given by 

T[3 F ] := iTrLn[S F ] - iTr[S^S F ] + T loop [S F }. (11) 

The last term of Eq. (|T1D is the residual term. Multiplying a factor i, iT\ oop [S F } is given by the 
sum of all Feynman amplitudes of 2-loop or higher- loop 2-particle irreducible vacuum diagrams 
in which every bare quark propagator 

S (x,y) = I e-*<'-»> - 1 - (12) 

h j{-q)9-™o 

is replaced by the full one 

S F (x,y) = (0\Tj>{x)$(y)\0). (13) 
The SD equation is the stability condition of the CJT action 

= 0. (14) 

SS Fmn (x,y) 

Throughout this paper, we employ the lowest order (lowest-loop) expansion of the ri oop [SV] as 

IWSf] = - l - j d A x1C m ^ n ^ (id Xl) id X2 -id yi) id y2 ) (15) 
x [S Fm2mi (x2, Xi)S Fn2Tll (t/2j Hi) ~ S Fm2ni (x2,yi)S Fn2mi (y2, X\)} 

where the symbol * means to taking x\,X2,yi,y2 —> x after all the derivatives are operated. 
This leads to the ILA model, where the SD equation is given in momentum space by 

zS^(q) - iS^iq) + C F [ f{qW E )iD^{jp - q)^S F (p) Jv = , (16) 

Jp 

with 

tr [T°T>] Nl-l 

Cf = ^^ = ^^- (17) 
Now we introduce the regularized propagator as 

S?(q) := f(-q 2 )S (q) = -r-^—, (18) 
S$(q) := f(-q 2 )S F (q). (19) 



Then the SD equation (|T6|) becomes 



l s*-\q)-<-\ q ) + 



c h 



?(q A E ,p E )iD^(p- ff) 7 ^#(p)7* = 



(20) 



/HzVp/H? 2 )' 

in which one finds that the integral is cut-off at p\ = —p 2 = Ayy due to the function f(—p 2 ). 
Substituting the general form of the SD solution 



A(q 2 )i-B(q 2 ) 
we obtain a set of integral equations 

iC F f 9 2 (q 2 E ,p 2 E ) 3(p 2 + q 2 )(pq) - A(pq) 2 - 2p 2 q 2 



(21) 



AW 



1 + 



x 



B{q 2 



m + 



q 2 f(-q 2 )Jp f(-p 2 ) 
A{p 2 ) 
p 2 A 2 (p 2 )- B 2 (p 2 Y 

iC F fg 2 (q 2 E ,p 2 E ) 1 



(q - pf 



-3B(p 2 



f(-Q 2 )Jp f(~P 2 ) (p-q) 2 p 2 A 2 (p 2 )-B*(p*y 
After the Wick rotation, we obtain 

A(-q 2 E ) = 1 



(22) 
(23) 



(24) 



from Eq. (|22"D . This is another advantage of the Higashijima-Miransky approximation, where 
the running coupling constant is defined so as to make the wave function renormalization Z 2 



unity j|, 10 1 . Then we find an integral equation for B(—q E ) as 

2, _ / / 7,2 -2/2 J\ PE B{~Pe) 



B (-qE) = m o + 



167T 2 Jo 



dp 2 E g 2 (q 2 E ,p 2 E ) 



max{g|, p%} p\ + B 2 (-p 2 E ) ' 



(25) 



4 Renormalization of quark mass 

In this section we discuss the renormalization of the quark mass. The operator product expan- 
sion analysis shows that in the asymptotic region the QCD quark mass function B(—q E ) for 
three quark flavors behaves as follows [O. 



9 2 (^ 



3q 2 



g_M) 



-4/9 



(26) 



where m^(/i 2 ) is the current quark mass renormalized at /i 2 and £,r{h 2 ) := {ip^R is the quark 
condensate renormalized at /i 2 . The improved ladder model of QCD is the model which is 



constructed so as to reproduce the QCD asymptotic behavior. Therefore we introduce the 
renormalization condition of the quark mass so that the solution of the SD equation can be 
interpreted as the QCD quark mass function in the asymptotic region. 

The quark mass function calculated in the effective model of QCD can be expressed in the 
similar fashion as Eq.(BBI) 

B{-q 2 E ) = m R ^ 2 ) F(q 2 E , p 2 ) - G(q%, p 2 ) . (27) 



Then we introduce the renormalization condition 



F(/i 2 ,/i 2 ^ 
which is equivalent to 

ae(-/i 2 ) 



dm R {jj?- 



1 . 



(2? 



(29) 



m R (n 2 )=0 



This mass independent renormalization condition for the SD equation is first proposed by 
Kusaka [11]. The mass renormalization constant Z m = Z m (Ay Y , /i 2 ) is introduced by 



™o( A uv) = Z m m R (fi ) 



(30) 



It should be noted here that the renormalization constant is independent of mass in this renor- 
malization scheme. It will be explicitly shown later in this section. By substituting itlq by 
Eq.flSOD, the SD equation (|25| ) becomes 

B(-p%) 



B(-q 2 E ) = Z-'mnifi 2 ) + f W dp 2 E K(q 2 E ,p 2 E ) 

Jo 



pI + b\- p 2 e 



2 \ ' 



(31) 



with 



K(q 2 E ,p%) := 



3Cf _2, 2 2\ 

9 {q E iPE) 



Pe 



167r 2- Vax{g 2 ,p|} 
By differentiating this equation with respect to m R (p 2 ) and taking m R (p 2 ) = 0, one obtains 
dB{-q%) 

m R (p?)=0 



(32) 



dm R (p?) 



Z^ 



+ 



A 2 
"uv 



MAM pl ~ Bl[ - pl) m ~ pl) 



(p 2 E + B 2 (-p 2 E )f dm R (j* 



.(33) 



m R (n 2 )=0 
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The renormalization condition fp9|) leads to 

(34) 



•A) (Pe + b o{-Pe)) omRK^i 



m R (ti 2 )=0 



Here f?o( — q E ) is the solution of the SD equation in the chiral limit, namely, 

2 , 3C F [ A UV 2—2/2 2 \ PE B (-p 2 E ) 



Bo(-q E ) = / dp%g\qiy E ) f^T 2 7r 2 7 2 V ^ 

i6vr 2 Jo max{g£,p|.} p| + B 2 {-p%) 



Now the combination of Eqs.(|33"D and (j3~4"D yields an integral equation for 9m ( g f j 

Eq.(p4|) explicitly shows that the mass renormalization constant Z m does not depend on the 



quark mass. 

In order to obtain the quark mass function with the renormalized current quark mass, one 
first calculates B (—q E ) by solving the SD equation in the chiral limit Eq. (|35|) . Next Z m and 
fff] are obtained by solving Eqs.([33|) and Eq. (|34]) . Finally Eq. (|3TD is solved to find 

B{-q E ). 

Let us now propose the following definition of the quark condensate. 

UV 2 ) ■■= Z-%(Aly) , (36) 
£o(A 2 UV ) := " ( / AUV l^rMS^q)] ~ f*™ i^rMS^qUt]) , (37) 



(2tt) 4 l ' v ^ j J (2tt 



= dS*(g) 
dm R (fi 2 ) 



m R {fi 2 ) . (39) 



In order to avoid the divergence originated by the perturbative quark mass contribution to 
the quark condensate in the UV region, the perturbative quark mass contribution should be 
subtracted. The key point of our definition of the quark condensate is Eq. (p9|) , namely, the 
perturbative quark mass contribution is defined using the fully calculated quark mass function 
B(q 2 ). The subtraction of the perturbative quark mass contribution obtained by the operator 
product expansion approach in the UV region is not sufficient. Our definition of the quark 
condensate has the desirable property: 

m (A 2 v ) £o(Auv) = m R (fi 2 ) ^(/i 2 ) . (40) 
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5 BS equation for Pseudoscalar Mesons 

The homogeneous BS equation is given by 

r q , ^Ti^ ( , , Xn> m '(y',x';P B ) = 0. (41) 
dS Fmn (x,y)dS F n'm'{y , x ) 

Here the BS amplitude is defined by 

Xnm(y,x-P B ) := (0\TiP n (y)^ m (x)\P B ) (42) 
for a q-q state \P B )- The normalization condition is (Pb\P'b) = (27t) 3 2Pbo^ 3 (-Pb — P'b) anc ^ 



Pb '■= (y M B + P b ,Pb) is the on-shell momentum. Eq.([H|) is expressed in momentum space 
S F \q + ) X (q; P B )S F 1 (q-) = -iC F f f(q 2 E , k%)iD^(q - k)^ X (k; P B )ju , (43) 

J k 

with 

, Pb (aa\ 

q+ = q + ~2, 1- = 1 - -j- , (44) 

where the Fourier transformation of the BS amplitude is defined by 

Xnm(y,x;P B ) = e~ iP * x [ e-^- x \ nm {q;P B ), X = *±±. (45) 



We introduce the regularized BS amplitude by 

XnmiV P) ■= f(-Ql)Xnm(q; P)f(-Q-), (46) 

then Eq.(^) is rewritten as 

sr\q + )x R (q;Ps)sr\q-) 

= ~ iCf I /( _ A .2 ) 1 /( _ A . 2) g 2 (gl. - kh,X R (k; P B ) lv . (47) 

We see again that the integral equation is regularized correctly. 

The BS amplitude for the pseudoscalar meson can be written in terms of four scalar ampli- 
tudes, 



v tf ( k . P \ _ x ( A °)g/ 



+ l (j)T (k;p)(rit-m)-r, 



s{k ] P) + <p P {k\P)k l + <P Q {k ] P)P 

{Ah 



bit 
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where X a denotes the flavor matrix. Substituting this into Eq. p7|) , we obtain coupled integral 
equations for four scalar amplitudes. The explicit form is rather complicated and given in 
appendix. The integral equations can be written down formally 

Mq; Pb) = I M AB {q, fc; P B )Mk; P B ) (49) 

J k 

where A,B denotes S,P,Q,T. Among the four dimensional integration d 4 k, two of the inte- 
grations can be performed analytically after the Wick rotation and we set the total momentum 
Pbe = (Me, 0, 0, 0). Then we obtain 



<Pa{qr, qe] M E ) = / dk R dk e M AB {q R , qe\ k R , k g ; M E )<j) B (k R , k e \ M E ) (50) 

J(k R ,k g )ei 

where 

k E = k 2 R , k E PE = k R M E sin k e (51) 
and the integral region is given by 
I := I (k R , kg) 



M 2 1 

k\ ± k R M E sin k e + — p < Auv \ ■ (52) 



4 

This integral region is determined uniquely by the cut-off function in Eq. Q4"7D and is consistent 
with the SD equation (p5|). 

There is no solution of Eq.([50"D for a real M E because M E is a Euclidean meson mass whose 
square is negative M E = —M B . Since the SD equation can be solved only for space like region 
q\ > 0, the region M E < is not accessible. Instead of solving Eq.(|50"l), we convert it into an 
eigenvalue equation for a fixed M E > 0, given by 

^<pA(qR, qe] M E ) = J dk R dk e M A B(qR, qe; k R , k d ; M E )(pB(kR, k e ; M E ) (53) 



where A is the eigenvalue that is equal to unity for the solutions of Eq.(|50|). We solve Eq.(^3|) 
numerically using the iteration procedure. When we iterate Eq. (|53"D , the eigen-f unction asso- 
ciated with the maximum absolute eigenvalue is dominated. Then we obtain the maximum 
eigenvalue and its eigen-function. 

In the numerical calculations, we use the discretization of the continuous variable (q R , qe) 
and come across a problem that the kernel K AB diverges at the point (q R , qe) = (k R , k 9 ). This 
divergence, which is originated from the gluon propagator, does not cause a real divergence. 
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We may remove this divergence by carefully choosing the discretization points in the iteration 
procedure. 

Once we obtain (the largest absolute) A as a function of M| > 0, then we extrapolate A(M|) 
to the time-like region M\ < and look for the on-shell point where A(— M§) = 1. Since this 
extrapolation is the most ambiguous procedure in our calculation, we will later consider another 
function which is similarly extrapolated to the one-shell value and compare the predicted values 
M\ obtained in the two independent extrapolations. 



6 Decay Constant and low-energy relation 

To obtain the decay constant, we need the normalization of the BS amplitude. The normaliza- 
tion condition of the BS amplitude is derived from the inhomogeneous BS equation 

^ , S ^! S Q F] , , . C w (^; x "v") = Sm»m5nn»5(x" - x)5(y - y") (54) 

where 

- (o\TMy)^ m (x)\o)(o\TMx')^M)\o)- (55) 

In the momentum space, Eq.(^) gives 

1 d 



-(^),/'(-g?.)/H£) - (<z- W (-g 2 )/(-£) 



+ ^ 



/ 2 (-^)/ 2 (-<? 2 ) 

Xn im M'-> PB)SF J ^ im2 {qJ)Xm2n 2 ( < l'-> P b) Sp n * ni (q + ) = -2P M , P — > P B . (56) 



In the case of the sharp cut-off function (|3|), the second term in the LHS of Eq.(^) does not 
contribute and the integral region in the first term is determined uniquely. 

Let us now turn to the discussion of the axial-vector Ward-Takahashi (WT) identity. It 
has been found that the axial-vector WT identity is violated in the Higashijima-Miransky 



(HM) approximation [TJJ]. Of course the Goldstone theorem holds in this case because the 



HM approximation respects the global chiral symmetry. The chiral WT identity in the ladder 



approximation has been carefully studied in ||10|| . The reason of the violation of the axial- vector 
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WT identity is that the HM approximation breaks local chiral symmetry. As shown in [fLC[| , 
the improved ladder approximation of the SD and BS equations preserves the WT identity for 
the axial- vector vertex if and only if one uses the gluon momentum square as the argument of 
the running coupling constant. However, in this case renormalization factor Z2 of the quark 
wave function deviates from unity in the Landau gauge. In order to avoid such problems, 



authors of [10] have introduced the non-local gauge so that the gauge parameter in the gluon 
propagator becomes a momentum dependent function. On the other hand, we have proposed 
another way to recover the axial- vector WT identity [JTJJ. In this approach the axial- vector 
current is modified so as to become the correct Noether current of the effective model of QCD. 
The advantage of this approach is that it is applicable to all the effective models of QCD 
which respect the global chiral symmetry. According to [T5| , the meson decay constant can be 
expressed as follows^. 



B 



lim — / tr 

P^P B P 2 Jq 



X a ff(-q 2 -) + f(-q 2 + ) 



2 V 
+E a (q;P) 



+ (f(-&) - f(-d) 



(57) 



- /C"'"< w (-k, q - |; -q + | , *) } (vysyM*) 

P P r 

-k + P,q-—;-q-—,k 



I ^ j^n n,mn 



-Kq + ^-q-^k)}h F (k)i l5 ^) 

\ / m'n' -1 



(5? 



The on-shell value /b(M| = — M|) is obtained again by extrapolation from the space like 
region M| > to the on-shell point M\ = —M\. For the neutral pion (B = tt°), a = 3 and so 
on in Eq. (E7f) . 

As shown in Ref . [|15|1 , the WT identity for the axial- vector vertex leads to the following 
relation in the improved ladder approximation model of QCD: 



M 2 B f B = -2m R ((x 2 )8 B (^ 



(59) 



/„. in this paper corresponds to f v in 15 
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with 



) := Z m (/i )z / ^tr 



2 

This relation is satisfied for a finite quark mass. In the case of chiral limit, we obtain 



(60) 



f B £ B (fx 2 ) = U^ 2 ) where £ B := (61) 

One can treat Eq. (|6lD as an approximated relation of the leading term of the expansion of m B for 
finite quark mass. Substituting Eq.(^) to Eq . fl59|) , one obtains the Gell-Mann-Oakes-Renner 
(GMOR) mass formula 

M lfl - -2m R {^) R | chiral Umit . (62) 



For finite m# > 0, Eq.(|59|) is an exact relation, while the violation of the GMOR formula is 
incurred by the violation of Eq. (|6T|) . 
We define 71 by 

, -Mif B (Ml) , . 

K {Ml) := f B \ EJ . 63 

-2m R £ B 

A relation TZ(— M|) = 1 must be satisfied due to Eq. (|59|) . We use this condition to make the 
extrapolation more reliable. 



7 Numerical results 

7.1 Parameters of the model 

The parameters of the improved ladder model of QCD are the current quark mass m B for up and 
down quarks (The isospin symmetry is assumed throughout this paper.), the scale parameter of 
QCD Aqcd, the infrared cut-off trp for the running coupling constant, the smoothness parameter 
to and the ultraviolet cut-off Auv- We take to = ~3 throughout this paper, which is the same 
value used in Ref. 0. In Ref.|Q they have shown that the numerical results are quite insensitive 
to the choice of the smoothness parameter in their studies of the BS amplitudes in the chiral 
limit. As for the ultraviolet cut-off Auv> we shall show that the physical observables depends on 
it rather weakly after our renormalization procedure described in section ^, if we use reasonably 
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large value of Auv- Of course, as we are treating not the full QCD but its approximation, we 
expect that weak dependences remain in our numerical studies. Thus we only have three 
physically relevant parameters, namely, the current quark mass, Aqcd and the infrared cut-off. 

We choose Aqcd about 0.6 [GeV]. It is rather large compared with the value obtained from 
the analyses in the deep inelastic scattering. In the framework of this model, however, one must 
employ the large value of Aqcd i n order to bring sufficiently strong dynamical chiral symmetry 
breaking. It may be the indication of the limitation of the improved ladder approach. Other 
non-perturbative interactions between quarks may solve this discrepancy. One candidate for 
such non-perturbative interactions is the instanton induced interaction proposed by 't Hooft 



T6|. There have been many studies of the roles of the instanton in low-energy QCD such as 



the instanton liquid model [I7|, the generalized Nambu-Jona-Lasinio (NJL) model |Tj| with 
't Hooft instanton induced interaction ||19|| , the effects of the instanton in baryon sector [pOfl . 
The recent studies of the //-meson properties in the generalized NJL model with the 't Hooft 
instanton induced interaction have shown [21] that the contribution from the 't Hooft instanton 
induced interaction to the dynamical mass of the up and down quark mass is 44% of that from 
the usual Ul(3) x Ur(3) invariant four-quark interaction. The introduction of the 't Hooft 
instanton induced interaction to the ILA model seems to be interesting and such attempt is 
now in progress p2| . 

We employ about —0.5 due to Ref.Q. The running coupling constant for various Aqcd 
and tip is shown in Figs.ffl and H. 



7.2 SD equation 

We discuss the solutions of the SD equation in this subsection. B(—q E ) as solutions of Eq. (f25|) 
for various values of m# are shown in FigFJ. In the chiral limit, i.e. ttlr = 0, we find a 
non-trivial solution B(—q E ) which is non-zero for q E < lGeV. Note that Eq.(p5"D has also a 
trivial solution B(—q E ) = 0, if = 0. The existence of the non-trivial solution indicates 
that chiral symmetry is broken dynamically. B(—q E ) decreases quickly to zero for qE > lGeV. 
The asymptotic behavior is consistent with the OPE result as shown in Ref.0. The range of 
non-zero value of B(—q E ) is determined by the q 2 E dependence of the coupling constant, g(q%) 
defined in Eq.(|7]). The results shown in Fig.|| correspond to the choice Aqcd = 0.6 GeV. When 
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Figure 1: q\ dependence of g 2 (q E ) for various A QCD with t IF = —0.5. 
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Figure 2: q 2 E dependence of g 2 {q 2 E ) for various t IF with Aq C d = 0.6GeV. 
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2.0 



Figure 3: Quark mass function B(—q E ) as function of qg. Model parameters are Auv = 100 
GeV, Aqcd = 0.6 GeV, t 1F = -0.5 and /i 2 = 4 GeV 2 . 





-0.5 0.0 0.5 1.0 1.5 2.0 


[MeV] 


259 240 187 116 43 2 



Table 1: tip dependence of quark condensate in the chiral limit. Other model parameters are 
Auv = 100 GeV, A QCD = 0.6 GeV and fi 2 = 4 GeV 2 . 

Aqcd decreases, the range of B(—q E ) decreases accordingly and therefore the chiral symmetry 
breaking is weakened. But the Aqcd determines the scale of the system. The order parameter 
(ipip) is almost proportional to the Aq CD and the chiral symmetry breaking always occurs for 
smaller Aqcd- The parameter tip determines the strength of the coupling constant. Table [I] 
shows tip dependence of the condensate in the chiral limit. As can be seen from table [1], the 
chiral symmetry breaking does not occur for larger tip. The asymptotic behavior of B(—q E ) 
with the finite current quark mass is rather different from that in the chiral limit. It can be 
seen clearly from the log-log plot of the quark mass function B(—q E ) given in Fig. |j. 

We next discuss the quark condensate. The quark condensates —{ipip/R calculated for the 
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Figure 4: Log-log plot of quark mass function B(—q%). Model parameters are A uv = 100 GeV, 
Aqcd = 0.6 GeV, t IF = -0.5 and /i 2 = 4 GeV 2 . 
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Figure 5: dependence of — (ipip) 1 ^ 3 . Model parameters are Auv 
GeV, % = -0.5 and /i 2 = 4 GeV 2 . 



100 GeV, A QCD = 0.6 



various quark masses are shown in Fig. |5]. Since in our definition of the quark condensate given 

1/3 

in Eqs.(|36|)-(|39|) the perturbative quark mass contribution is subtracted, —{ipip)R decreases as 
uir increases. Similar behavior is observed in the QCD sum rule approach. {ip^R for ttlr = 120 
MeV is about 78% of (ipip)n for = 5 MeV. This is in reasonable agreement with the QCD 



sum rule result p3|: (ss)/{uu) = 0.8 ± 0.1. 



7.3 BS equation 

Let us now turn to the discussion of the solutions of the BS equation. The eigenvalues A(M|) 
of the BS equation are shown in Fig.|| One sees that the massless solution A(Mg = 0) = 1 
appears in the chiral limit. This is a result of the Nambu-Goldstone theorem. For a non-zero 
quark mass, we need to extrapolate A(M|) to the time-like M§ < 0. We have fitted A(M|) by 
a quadratic function using the method of least-squares and extrapolated A(M§) to the time-like 
region to find the point at which A(Mj,) becomes unity. Fig.^ implies that the extrapolation 
length is longer for larger quark mass. In order to reduce the ambiguity in the extrapolation 
procedure, we also evaluate the ratio 71 defined by Eq.(^) as a function of M|. Because of the 
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I m R =0 MeV 

I m R =5 MeV 

I — — m R =10MeV 

I m R =20 MeV 

I 
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-0.10 0.00 0.10 0.20 

M| [GeV 2 ] 

Figure 6: Eigenvalues of BS equation A(Mf). Model parameters are Auv = 100 GeV, Aqcd = 
0.6 GeV, t w = -0.5 and /i 2 = 4 GeV 2 . 

exact relation (^), 7Z must hit TZ(-M^) = 1 at the on-mass-shell of the pion. The value TZ is 
shown in Fig.[7|. We fit the graph with the linear function using the methods of least-squares. 
We show our calculated results of the pion mass M n determined by the above mentioned two 
conditions, i.e., A = 1 and TZ = 1, for various values of the quark mass m R in table |2|. The 
ambiguity by the extrapolation procedure is reasonably small up to the strange quark mass 
region. We plot our results of M 2 obtained by the condition A = 1 as a function of m R in Fig. 
||. The M 2 seems to be almost a linear function of m R up to m R ~ 40 MeV. This is suggested 



m R 


5 10 20 40 80 120 


M n (A = 1) 


0.0 159.1 222.0 312.9 444.4 639.8 800.9 


M n (TZ = 1) 


0.0 154.5 218.5 309.0 436.9 616.1 749.7 



Table 2: Pion masses determined by the two conditions, A = 1 and TZ — 1 for various values 
of rriR. All the entries are in units of MeV. Model parameters are A uv = 100 GeV, Aq C d = 0.6 
GeV, % = -0.5 and /i 2 = 4 GeV 2 . 
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Figure 7: The value TZ(M^). Model parameters are Ayv = 100 GeV, A QCD = 0.6 GeV, 
t lF = -0.5 and /i 2 = 4 GeV 2 . 

by the GMOR formula 

J2 J . ( 64 ) 

J n / chiral limit 

The deviation from the linear dependence at = 120 MeV is about 9%. 

We next discuss the pion decay constant. As mentioned in Sec. [| we can calculate /^(Mf) 
only for the time-like M\ and therefore the on-shell value of the decay constant /^(Mf = — M 2 ) 
can be obtained again by the extrapolation. The M\ dependence of the pion decay constants 
in the chiral limit and in the case of m# = 5 MeV are shown in Fig. [J[ 

To estimate the effect of E a (q; P) in Eq . (|57|) , we plot the naive value which is defined 

by neglecting E a (q; P) term from Eq.([57|). It seems to be a good approximation that f n (M^) 

is a linear function of Mf. Therefore we fit the curve by the linear function using the method 

of least-squares and make an extrapolation to the time-like Mj, for finite m Q. On the other 

hand, we fit by the quadratic function using the least-square method to extrapolate to the 

It should be noted that the decay constant at small (positive) Mjj suffers from numerical uncertainty and 
thus it deviates from the straight line. We do not use these points in our extrapolation procedure. 
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Figure 8: M 2 as a function of m R . Dots represent our calculated results and line represents 
the linear fit of the first 4 points. Model parameters are A uv = 100 GeV, Aq CD = 0.6 GeV, 
t w = -0.5 and /z 2 = 4 GeV 2 . 
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Figure 9: M\ dependence of f n and f^. Vertical dashed line represents the point M| = — M 2 
for = 5 MeV. Model parameters are Auv = 100 GeV, Aqcd = 0.6 GeV, tip = —0.5 and 
^2 = 4 GeV 2 . 

on-shell point M| = M 2 . Our results for m R = 5 MeV are f v = 88 MeV and f£ = 155 MeV, 
so the contribution of E a (q; P) term is remarkable. Similar result has been found in the ILA 
model which respects the axial-vector WT identity by using the gluon momentum square as 



the argument of the running coupling constant and the non-local gauge flO|| . The condition 
K{-Ml) = 1 is the direct consequence of the axial- vector WT identity and therefore it has been 
proved numerically that our definition of the decay constant is consistent with the axial-vector 
WT identity from the fact that the pion mass determined by the condition A = 1 is almost 
same as that determined by the condition 1Z — 1. 

We plot the quark mass dependence of the decay constant in Fig. [Tl]. /„. almost linearly 
depends on m R in our case. f n at m R = 40 MeV is about 15% bigger than f n at m R = 5 MeV 
and f n at m R = 80 MeV is about 32% bigger than f n at m R = 5 MeV. Since the observed 
fx/ fn = 1-23, m R dependence of the decay constant seems to be reasonable though we have 
not solved the kaon BS equation. This m R dependence of /„- is similar to that obtained in 
the chiral perturbation theory (ChPT) |24| though the chiral log term has not been seen in 



our numerical result. It is understandable because the Goldstone boson loop contribution is 
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Figure 10: m R dependence of /„.. Model parameters are Auv = 100 GeV, A QCD = 0.6 GeV, 
tjp = -0.5 and [i 2 = 4 GeV 2 . 



Aqcd [GeV] 


M w [MeV] 


U [MeV] 


-(V^)i /3 [MeV] 


0.5 


152 


74 


223 


0.6 


159 


88 


259 


0.7 


166 


102 


293 



1/3 

Table 3: A QCD dependences of M n , f n and —{4>?P)r ■ Other model parameters are Auv — 100 
GeV, t 1F = -0.5, [i 2 = 4 GeV 2 and m R = 5 MeV. 

not taken into account explicitly in our approach. On the other hand, the quark-antiquark 
structure is included explicitly and the finite quark mass effects are fully taken into account 
without performing the perturbative expansion with respect to the quark mass. 

Let us now discuss the Aqcd and the tw dependences. Table ^| shows the Aqcd dependences 
of the pion mass, the pion decay constant and the quark condensate. As shown in Table [3], 
the all the quantities with the mass dimension one are roughly proportional to the Aqcd- It 
is understandable since the only scale of the theory is the Aqcd if one can neglect the current 
quark mass. Table f| shows the tw dependence. For below —0.7 the coupling constant 
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tlF 


M n [MeV] 


U [MeV] 


[MeV] 


-0.3 


149 


91 


256 


-0.5 


159 


88 


259 


-0.7 


(181) 


(74) 


(253) 



1/3 

Table 4: i IF dependences of M^, /„. and — {ipip) R ■ Other model parameters are A uv = 100 
GeV, Aqcd = 0.6 GeV, /i 2 = 4 GeV 2 and m R = 5 MeV. 

becomes very steep and our numerical procedure is not sufficiently accurate. Although we have 
not performed the fine tuning of the model parameters, it is clear from Tables [l], ||] and £| that 
one can fit the model parameters so as to reproduce the observed values of M n and f n and the 
empirically determined value of the quark condensate. 

We study the Auv dependence by changing the value of Auv from 10 GeV to 1000 GeV. 
It causes less than 1% changes of the M^, f n and (ipip)ji. This stability indicates that our 
non-perturbative renormalization procedure works well. 



7.4 Approximation 

Finally we discuss the approximation often used in solving the BS equation. The approximation 



in which one neglects the 0p,0q,0x terms in RHS in the BS equation ((49|) is often used in 



literatures p5| . The resulting eigenvalues are shown in Fig.|ll|. While this approximation gives 
the massless NG boson in the chiral limit, it underestimates the pion mass for finite quark 
mass. The decay constants obtained from the approximated BS amplitude are shown in Table 
H The approximation overestimates the pion decay constant about 30% while the axial-vector 
WT identity (M^ ppro -) 2 / Appro - = -2m R £^ ppro - is preserved. It is seen that the effect of E a (q; P) 
is very small in this case, and f^ ppTO ' — f^ A PP TO -. Therefore, the violation of the axial- vector 
WT identity or that of the exact relation for the PCAC current incured by neglecting E a (q; P) 
effect is very small. 

The above f^ ppro ' is calculated in a similar way as in the approximation discussed in Ref. Hi~5 



26| . The following is shown in theorem 2 of [pqj . If the interaction is local chiral invariant, 
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Figure 11: Eigenvalues of the BS equations with and without the approximation described in 
the text. Model parameters are Auv = 100 GeV, Aqcd = 0.6 GeV, tpF = —0.5 and /z 2 = 4 
GeV 2 . 



m R 


M„ 


nf Appro. 

7T 


u 


yAppro. 


fN 


yiVAppro. 











86 


115 


154 


116 


5 


159 


119 


88 


117 


156 


117 


10 


222 


166 


90 


118 


158 


119 



Table 5: Result with and without Approximation. All the entries are in units of MeV. Model 
parameters are A uv = 100 GeV, A QCD = 0.6 GeV, tip = —0.5 and /i 2 = 4 GeV 2 . 
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the approximation of taking the wave function renormalization of the quark propagator to be 
one and at the same time neglecting the (ftp, <pQ and (ftx terms in the RHS of the BS equation 
preserves the low-energy relations. This theorem cannot be applied to the present case since 
the interaction term of the ILA model breaks the local chiral symmetry. However one can prove 
that the low-energy relation holds if one neglects the (ftp and <pQ terms in E a (q; P) as well as 
the (ftp, (ftn and (ftx terms in the RHS of the BS equation by following the same argument of 



the proof of theorem 2 in [26 



8 Conclusion 

We have solved the Schwinger-Dyson (SD) equation for the quark propagator and the Bethe- 
Salpeter (BS) equation for the pion in the improved ladder approximation of QCD. We have 
carefully treated the consistency of the equations in order to preserve the low-energy relations 
associated with chiral symmetry by using the Cornwall- Jackiw-Tomboulis effective action ap- 
proach. We have introduced the finite quark mass term in order to study effects of explicit chiral 
symmetry breaking on the low-energy relations. Because of the difference in the asymptotic be- 
havior of the quark mass function for finite quark mass, the non-perturbative mass-independent 
renormalization has been introduced and the quark condensate for finite quark mass is calcu- 
lated. In solving the SD and BS equations, we have not taken any further approximation such 
as expansion of BS amplitudes in the Gegenbauer polynomials. 

We have obtained reasonable values of M n , f n and {i/)iP)r with a rather large value of Aqcd- 
It may indicate the limitation of the improved ladder approach. The pion mass grows as 
quark mass mp increases. Up to the strange quark mass region M% seems to be proportional 
to quark mass mp almost as predicted by the GMOR relation 

M * = 75 m R- ( 65 ) 



j - 



chiral limit 



We have found that the f n also grows as mp increases almost linearly. The mp dependences of 
M% and f n are similar to those obtained in the chiral perturbation theory. It suggests that the 
chiral perturbation is applicable up to the strange quark mass region. 

We have studied the effect of E a (q;P) term in the true decay constant. We have found 
that it is significantly large for various input parameters. Therefore in the framework of the 
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improved ladder approximation, E a (q; P) plays an essential role to keep the chiral property. 

We have further shown the result of the approximation neglecting (fip(q; P), (f)Q(q; P) and 
4>t(q', P) term in RHS of the BS equations. This approximation is very useful and makes the 
calculation easy greatly. But the result gives a smaller pion mass. This suggests that the simple 
picture of the (ps(q',P) dominance in the BS equation is not so good, at least in the present 
model. 

So far, we have studied the symmetric q-q systems, uu, dd, etc. It is interesting to extend 
the present formulation to asymmetric systems like the kaon. It is also interesting to introduce 
the Ua(1) breaking interaction to this framework and to study the r\-r( systems. Such attempts 
are in progress. 
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Appendix 

Here we write down the BS equation explicitly. In this section the total momentum is denoted 
by P instead by Pb for simplicity. First we define the regularized amputated BS amplitude 



X R (q; P) by 

X R (q; P) := Sr\q + ^)x R (q; P)S§r\q ~ \ ) 



(66) 



which can be expressed in terms of 




(67) 



The BS equation (|49"D reads 




(68) 
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The components of the kernel is given explicitly by 

-3 



Kss(q,k;P) 
K PP (q,k;P) 



iC F g (q,k)- {q _ k)2 
iC F g 2 (q,k) r P 2 (qk) - (Pq)(Pk) 



(69) 



K PQ (q,k;P) 
K QP (q,k;P) 



K QQ (q,k;P) 



P 2 q 2 _ ( Pg )2 ( g _ k y 

2{qk - k 2 )(P 2 q 2 - {Pq) 2 + {Pq){Pk) - P 2 {qk)) 1 
+ {q-kY J 

iC F g 2 (q, k) 2{Pq - Pk)(P 2 q 2 - (Pq) 2 + (Pq)(Pk) - P 2 (qk)) 
P 2 q 2 _ ( P? )2 (g - A;) 4 

tC F g 2 (q,k) f (Pk)q 2 - (Pq)(qk) 
P 2 g 2 -(Pg) 2 l (q-k) 2 

2{qk-k 2 ){{Pq){qk)-{Pk)q 2 ) \ 
(q- k) 4 J 
iC F g 2 (q,k) (P 2 q 2 -(Pq) 2 



(70) 
(71) 



(72) 



PV - (Pq) 2 I (g - A;) 2 

2(Pq - Pk)((Pq)(Pk) - (qk)q 2 ) j 



+ 



K TT (q,k;P) = 



(q - k) 4 
iC F g 2 {q,k) 1 



(73) 



(A; 2 -g 2 )((Pg)(P£;)-P 2 (#)) 



P2 g2 _ ( P? )2 ( 5 _ k y 
+ 2(Pg - Pk)((Pk)q 2 - (Pq)(qk)) - 2(qk - k 2 )(P 2 q 2 - (Pg) 2 )| (74) 



and other components are zero. The relations between 0^(g; P) and 0^(g; P) are given by 



Mq;P) = 



i 

A 



P 2 



5 



P(^)P( g !)U 5 (?;P) 



+ {g 2 (P(g 2 ) - B(q 2 _ )) - ^(B(q 2 _) + B(q 2 + ))} P (g; P) 

+ " ^(g 2 )) - ^{B{q 2 _ ) + P(g 2 ))U Q (g; P) 



+ (p 2 q 2 -(p q ) 2 )Mq;P) 



(75) 



(B(q 2 + )-B(q 2 _))4> s (q;P) 

f p2 ^ 



+ 2(Pq)4> Q (q;P) 



>(q;P) 



- (Pq)(B(qi) + B(q 2 + )) + -y (B(<£) - M<ZI P) 



(76) 



l 

A 



--(B(q 2 + )+B(q 2 _))Mq;P) 
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-^*p(r,P) 

-jg 2 + ^+P(^)P(g 2 )}<Mg;P) 

Pq 



+ {q 2 (B(ql) + B(qD) + -£(B(<?_) - B(<? + ))} ^(9; P) 



(77) 



+ (P(gi) + P(d))0,(g;P) 



+ \~ q2 + T ~ B(q 2 _)B(ql)\ T (g; P) 



(78) 



where 

■>2 / _,2 \ \ / Jl 7^)2 / 2 



A:=«-SV))(9l-^(9l))- (79) 
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